Fractional topological insulators of Cooper pairs induced by proximity effect 
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Certain insulating materials with strong spin-orbit coupling can conduct currents along their 
edges or surfaces. This phenomenon arises from the non-trivial topological properties of the elec- 
tronic band-structure, and is somewhat similar to the integer quantum Hall effect of electrons in 
strong magnetic fields. Topological insulators analogous to the fractional quantum Hall effect are 
also possible, but have not yet been observed in any material. Here we show that a quantum well 
made from a topological band insulator such as Bi2Se3 or Bi2Te3, placed in contact with a super- 
conductor, can be used to realize a two-dimensional topological state with macroscopic many-body 
quantum entanglement whose excitations carry fractional amounts of electron's charge and spin. 
This fractional topological insulator is a "pseudogap" state of induced spinful p-wave Cooper pairs, 
fundamentally significant to quantum field theory and applicable to spintronic devices and quantum 
computing. 



The recently discovered two-dimensional topological 
insulators (TI) with time-reversal (TR) symmetry [l|- 
Q] are band-insulators related to integer quantum Hall 
states in which electron spin plays the role of charge. 
They can be obtained in HgTe, Bi2Te3 and Bi2Se3 quan- 
tum wells owing to the strong spin-orbit coupling, and 
exhibit topologically protected gapless edge states de- 
spite the spin non-conservation [fj. The properties of 
quantum wells are linked to the topologically protected 
surface states of the extensively studied bulk materials 



Instabilities caused by interactions among electrons 
can establish unconventional quantum states in TIs, with 
broken symmetries 0, [l(| or topological order 11 - 14 1 . 
These envisioned forms of quantum matter could realize 
robust macroscopic entanglement between spatially sepa- 
rated electrons in the TI materials, which motivates both 
the fundamental research and the quest for applications 
in spintronics and quantum computing. Here we aim to 
realize a new class of strongly correlated TIs that ex- 
hibit phenomena reminiscent of the fractional quantum 
Hall effect (FQHE) in strong magnetic fields, but with- 
out its TR symmetry violation 1^-22|. Such fractional 
TIs feature quasiparticles that carry fractional amounts 
of electron's charge and spin. Exotic states with non- 
Abelian statistics are also possible and promise the abil- 
ity to perform quantum computation with a greater level 
of quantum control than in FQHE qubits, because both 
charge and spin can be manipulated and entangled. 

One approach to obtaining fractional TIs, inspired by 
the FQHE, exploits Coulomb interactions among elec- 
trons in a partially populated band made narrow by the 
spin-orbit coupling [23l425j . It might be very difficult 
to find TI materials with sufficiently narrow bands and 
strong interactions, so the goal of this paper is to propose 
a different approach. Here we consider a heterostruc- 
ture device in which a two-dimensional electron gas can 



be tuned near a quantum critical point. Every quan- 
tum critical system is sensitive to relevant perturbations 
that impose their energy scales on dynamics and define 
the quantum phases that surround the critical point in 
the phase diagram. We will show that the spin-orbit 
coupling is characterized by a large "cyclotron" energy, 
and thus indeed represents a relevant perturbation that 
can dominate near the critical point and stabilize frac- 
tional topological states just like a strong magnetic field 
would. The proposed heterostructure is not only rou- 
tinely achievable, but also provides the best platform to 
experimentally seek a variety of topologically non-trivial 
superconducting and insulating quantum states that have 
not been observed or hypothesized before, and whose ex- 
istence is guaranteed by the fundamental principles dis- 
cussed here. 

We engineer a quantum critical point by placing a 
TI quantum well in contact with a conventional super- 
conductor (SC) as shown in FigO] The SC's phonon- 
mediated pairing glue and electron tunneling into the 
TI generate the proximity effect, a short-range attrac- 
tive interaction between electrons inside the quantum 
well. This produces a superconducting state in the TI 
when its energy gap is reduced somewhat below the SC's 
critical temperature T c by applying a gate voltage. The 
critical gap is actually finite because the short-range at- 
traction helped by the TI's two-dimensionality is strong 
enough to create bound-state Cooper pairs. Therefore, 
the superconductor-insulator transition in the well, tuned 
via the gate, is driven by soft collective modes of Cooper 
pairs whose energy lies below the pairing gap. Some pairs 
are spinful triplets allowed by the TIs orbital degrees of 
freedom. Since the spin-orbit coupling significantly af- 
fects their dynamics near the transition, the existence of 
fractional topological states is assured in the "vortex liq- 
uid" region of the phase diagram in FigJ^a). Our goal 
here is to paint the physical picture of this using funda- 
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FIG. 1: The heterostructure device that can host fractional 
TR-invariant quantum states. A topological insulator (TI) 
quantum well is sandwiched between a conventional super- 
conductor (SC) and a conventional insulator (I). The SC's 
pairing glue and the virtual tunneling of Cooper pairs from 
the SC to the TI can induce two-dimensional superconduc- 
tivity or a strongly correlated topological insulator of Cooper 
pairs in the TI. The gate (G) voltage can be used to control 
the state of the TI, and the topological properties of the TI 
can be probed via a Hall-bar setup of leads (L). 

mental principles rather than a microscopic calculation. 

The minimal model Hamiltonian of a non-interacting 
TI quantum well is: 

H = v z(S x p)r z + At x - n . (1) 

It describes four electron states per momentum p, la- 
beled by the spin projection S z = ±| and the orbital 
index r z = ±1. We use the units with H = 1; the vec- 
tor spin operator is S = ^a a r a , a £ {x,y,z}, and a a 
and T a are Pauli matrices that operate on the spin and 
orbital states respectively. In this representation of the 
Bernevig, et al. model [2j, the orbital index is equivalent 
to the top/bottom surface of the quantum well. A ^ 
implements inter-surface tunneling that opens a gap in 
an otherwise massless Dirac spectrum of the spin-orbit 
coupling v. The model has relativistic particle-hole sym- 
metry when fi = 0. The many-body ground state of ([1]) 
is a band-insulator for \fj,\ < |A|, which is either trivial 
or topological depending on the sign of A and the high- 
energy spectrum that is not captured by ([I]). 

In order to expose the topological properties of this 
model we will compare it to a gauge theory of electrons 
in a static external Yang-Mills gauge field A [27] : 

H g = (P ~ TZA)2 + At x - n' = H a + ^- + const . (2) 
2m 2m 

The SU(2) gauge field A = —mv(z x S) embodies the 
spin-orbit coupling and has a non-zero "magnetic" flux 
density despite being uniform: 

= e^iduAx - ir z A v A\) = ^(mv) 2 S» T z a* . (3) 



(a) (c) 




FIG. 2: (a) The qualitative zero-temperature phase dia- 
gram of attractively interacting electrons in a two-dimensional 
band-insulating quantum well. Gate voltage Vg controls the 
electron gap in the quantum well and tunes the quantum 
phase transition (filled dot) between a superconductor (SC) 
and a bosonic Mott insulator (MI) of Cooper pairs. The low- 
est energy excitations in the MI are charge 2e bosons, but 
they become incoherent at the crossover (dashed line) to the 
band-insulator (BI) where only gapped fermionic quasiparti- 
cles with charge e exist. A spin-orbit coupling whose strength 
is measured by a "cyclotron" energy ui<s, (defined in the text) 
introduces a vortex lattice in the superconducting state of 
spinful p-wave Cooper pairs. Quantum melting of such a 
vortex lattice gives rise to correlated "vortex liquid" (VL) 
states, which are the prime candidates for fractional TIs. (b) 
The energy spectrum E(k) of the non-interacting SU(2) gauge 
Hamiltonian © for m ~ 2.5 • 10" 31 kg, v ~ 4 • 10 5 m/s and 
A ~ 100 eV (k — p/h). This two-orbital example approxi- 
mates the ARPES spectrum from the Fig.3d of Ref. |2q] when 
A = 0. (c) The Cooper pairing channels include intra-orbital 
spin singlets ($±), inter-orbital spin singlets (4>o), and inter- 
orbital p-wave spin triplets (rj m , where m £ {0, ±1} is the 
spin projection quantum number). 

Gauge transformations A^ — > WA t _ l W^ l +iT z Wd t _ l W^ and 
W^W^ are carried out by SU(2) matrices W = 
exp (|cr a 6' a ) , where 6 a , a 6 {x, y, z} are three coordinate- 
dependent angles. The "magnetic" flux $ = (mv) 2 and 
the corresponding cyclotron energy u)<& = mv 2 are the 
source of all quantum-Hall-like phenomena in TIs. 

The models {1} and ([2]) produce exactly the same dy- 
namics for single-particle charge currents j+ = i[H,Xi], 
as revealed by the operator equation of motion in the 
Heisenberg picture (at t = 0): 

^ = i[H, ji] = \v 2 e ijPj a z - Ave^r* . (4) 

The first term on the right-hand side reveals that clas- 
sical electrons with spin S = S z z would move along cy- 
clotron orbits whose sense of rotation (chirality) depends 
on their spin orientation. This fundamental source of all 
TI's topological properties is faithfully captured by H g , 
thus enabling the gauge theory description. Note that 
chirality is r 2 -independent, and hence unaffected by a 
finite bandgap A that only mixes the orbital content. 

The spectrum of @ is different from the massive Dirac 
spectrum of |T|) only by the Dirac cone curvature at 
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larger momenta, due to the extra p 2 /2m term. A small 
curvature is indeed seen in ARPES measurements [26j | . 
FigHJb) shows that §Z§ adequately approximates mate- 
rials, with a relatively large fitted mass m. The gauge 
symmetry of @ is nothing but a useful mathematical 
tool built on the approximate global spin SU(2) symme- 
try at low energies that is reduced to U(l) only by the 
spin-orbit coupling (<£> ^ 0). It will guide us through the 
analysis of Cooper pair dynamics in the TI, but at the 
end we will need to consider the consequences of gauge 
symmetry violations due to the spin non-conservation. 

The electron dynamics in the TI quantum well is al- 
tered by the proximity to the SC in the device from 
FigfT] The SC is a fully gapped quantum liquid of 
Cooper pairs characterized by two energy scales, the pair- 
ing A p and photon A 7 = hcX^ 1 gaps, where c is the 
speed of light and Al is the London penetration depth. 
Fermionic quasiparticles have anomalously small or van- 
ishing density of states below the pairing gap, which 
is A p = 2hvf/iT^ in conventional superconductors with 
Fermi velocity Vf and coherence length £. The smaller 
of the two gaps defines a cut-off energy scale for the low- 
energy dynamics in the TI that we shall discuss. The 
dynamics responsible for the superconductor-insulator 
transition in the TI is indeed defined below this cut-off 
and hence can be captured by a two-dimensional effec- 
tive theory whose degrees of freedom are decoupled from 
those of the SC. The fundamental U(l) symmetry allows 
proximity-induced attractive interactions in the effective 
theory, but prohibits explicit mean-field pairing terms. 
The SC's energy gaps then define a finite length-scale 
for the short-range attraction among the TPs electrons. 
We will show later that the resulting theory features a 
superconductor-insulator transition inside the TI across 



which A p ^ 0, A n 



0. 



Our effective TI model is given by the action S = 
J dt dr 2 (do + H g )ip + Si nt of the field-theory path inte- 
gral z = J "Dip e~ s defined in imaginary time t. The most 
important proximity-induced attractive interactions are: 



S\ n t — 



dtd 2 r 



(t/x^^VwVv. (5) 



where r = ±1 and a — ±1 label the orbital 
and S z states of the electron fields ijj TIJ respectively 
(r = — t). The unscreened Coulomb and orbital- non- 
conserving interactions are denoted by dots. By applying 
the Hubbard-Stratonovich transformation exp (§ie 2 ) oc 

J dy exp ^— |^ + xyj on the path-integral, we can elimi- 
nate the interaction couplings ([5]) in favor of six Cooper 
pair fields viewed in FiglJJc): two intra-orbital singlets 
4>± (Ui), two inter-orbital S z = ±1 triplets rj± (U 2 /3 at 
a=o J ), inter-orbital singlet <fio and S z — triplet rjo (U 2 /3 





t/v^k) 


<Mk) 


?7m(k) 


% translations 


Vv CT (k) 


<Mk) 


??m(k) 


IZg rotations 


VVcr^ek) 




r/ m (llgk.) 


IZi spatial reflect. 


V>Ttr(%k) 






Xt time reversal 


<W4(-k) 




(-l) m ^(-k) 


C charge U(l) 


e"V TCT (k) 


e l2 V„(k) 


e l2 % m (k) 


S spin U(l) 


e lCT >™(k) 


<Mk) 


e i2me Vm (k) 


local spin SU(2) 




<Mk) 


U mm' Vm' (^0 



TABLE I: Symmetry transformations of electron f/v T , singlet 
4> n and triplet r\m fields in (|6). The local SU(2) transforma- 
tion is relevant for the gauge theories ((2)1 and ([7]). W and 
U are SU(2) transformation matrices exp (i^/ a 9 a ) with SU(2) 
generators j a , a e {x, y, z} in the S = | and S = 1 represen- 
tations respectively, (r, a = ±1; n, m £ {±1, 0}; I = —I) 
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( U \<t>r\ 2 + 4> r Z„„4\A\a, + h.C 

=±1 

<T=±X 

+ y'\v \ 2 + v ^= (vVA + ^V-t) + hx 



The symmetry transformations of these fields are sum- 
marized in the Table [I] 

The condensation of any Cooper pair field in © 
corresponds to a superconductor-insulator quantum 
phase transition. Fermionic excitations remain gapped 
across such transitions in simple two-dimensional band- 
insulators with attractive interactions [28j j . and we will 
explain shortly why this also holds in TIs. Then, we 
may integrate-out the gapped fermion fields in the path- 
integral: 



to obtain a purely bosonic effective action S e g that de- 
scribes Cooper pair dynamics at energies below the pair- 
ing gap. Here we can rely only on symmetries to con- 
struct the Landau- Ginzburg form of S e s. Since two elec- 
trons with the same spin from different orbitals have the 
same cyclotron chirality, the Cooper pairs r/± with S z = 
±1 possess the SU(2) charge, unlike the singlet fields. 

Together with T70 they form a triplet r\ = (77 , 770 , ?7-i- ) 

that minimally couples to the spin-orbit SU(2) gauge field 
A. = —mv(z x S) expressed in the 5 = 1 representation. 
This can be seen from the local SU(2) transformations 
shown in the Table |TJ If all Cooper fields have modes at 
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energies below the two-electron continuum, we obtain: 



5 e g = dtdr 



(7) 



+U. 



S,(Ji(T2(T30-4 

The dots denote Coulomb interactions, and we omit- 
ted the non-relativistic corrections due to m < oo and 
H 7^ 0. We also organized the singlet fields into a vec- 
tor <p — ((/>_, 0o, (f>+) and expressed their quadratic cou- 
plings in the matrix form. The tensors K s ,t s ,t' s ,U s are 
only required to obey the TR symmetry, but they are 
not important for further discussion. The global SU(2) 
spin symmetry can be reduced only via the SU(2) flux in 
the gauge theory, while more realistic models allow other 
triplet field couplings. 

In a weak interface between a conventional SC and the 
TI, the SC's phonon mechanism is the main source of 
the induced attractive interactions among the TI's elec- 
trons. It originates in spin-insensitive Coulomb forces 
whose screened range is comparable to the quantum well 
thickness. Therefore, inter-orbital triplet pairing is pos- 
sible in the TI, but competes with singlet pairing. One of 
the intra-orbital singlet channels has a stronger induced 
interaction than all inter-orbital channels for geometric 
reasons, which naively means that singlets should con- 
dense before triplets when electrons are drawn into the TI 
from the SC by the gate voltage. However, the spin-orbit 
coupling mixes the spinful triplets into two helical modes, 
analogous to the Dirac particle- and hole-like eigenstates 
of @, one of which has energy that decreases when its 
momentum grows. That helical mode is destined to "al- 
ways" condense at a sufficiently large momentum accord- 
ing to , and thus has a natural advantage over singlets 
despite its origin in the weaker induced interaction. 

Being shaped by the spin-orbit SU(2) flux that can- 
not be expelled, a helical mode condensate is a peculiar 
yet to be observed superconductor with a lattice of lo- 
calized SU(2) topological defects. If the TR symmetry 
is not spontaneously broken, this is an Abrikosov lattice 
of coinciding quantized t] + vortices and r\~ antivortices 
without any spin texture, shown in Fig|3] Its existence 
also gives birth to fractional TIs. Imagine tuning any 
parameter that reduces the superconducting stiffness p s 
toward zero. The vortex kinetic energy due to zero-point 
quantum motion can be estimated from the Heisenberg 
uncertainty as -Ekin ~ ti v j where /$ is the SU(2) 
"magnetic length", and m v is the effective vortex mass. 
In (charged) superconductors, m v is roughly constant as 
p s — » 0, but turns into m v ~ | log(p s )| when the screen- 
ing length Al ~ ps 1 ^ diverges The potential energy 
due to the vortex lattice stiffness scales as E pot ~ p s per 
vortex (p 2 if the spectrum has Landau levels) , which is a 
simple estimate from the free energy expansion in pow- 




FIG. 3: A TR-invariant vortex lattice. Coinciding 77+ vor- 
tices and n~ antivortices are induced by the SU(2) spin-orbit 
flux. The counter-circulating supercurrents of 77+ and r\- are 
shown in red and blue respectively. Quantum melting of such 
a vortex lattice yields a fractional TI. 



ers of p s [30(. There is a critical finite p s at which the 
vortex lattice melts in a first order transition because 
£-kin > E pot . This happens at the solid line that sepa- 
rates the SC and VL regions in Fig|3Ja). Since p s also 
measures the quasiparticle pairing gap, its finite value 
implies that the transition is shaped by the Cooper pair 
dynamics below the fermion excitation scales. The result- 
ing insulator is a quantum vortex liquid of uncondensed 
Cooper pairs, whose qualitative properties are captured 
by the purely bosonic theory ([7|). We stress here that 
the presence of this regime is guaranteed by the general 
physical principles outlined in this paper. 

Quantum liquids of SU(2) topological defects are the 
prime candidates for fractional TI states when their den- 
sity is larger than the density of mobile particles. The 
mass m in ([2]) can be estimated from the curvatures of the 
top and bottom Dirac cones in ARPES experiments (2(| , 
and it is larger than the "spin-orbit" mass m so — A/v 2 
by a factor of A = ^- « 5 - 10 (v « 5 • 10 5 m/s). The 
cyclotron energy w$ = $/m = A A is not small in quan- 
tum wells with bandgaps A = 10 — 100 meV that can 
be engineered with a few atomic layers The den- 
sity of "magnetic" SU(2) flux quanta is rt$ = $//i 2 = 
X 2 A 2 /(vh) 2 w A 2 x 2 • 10 15 m" 2 . These estimates look 
promising if we compare them with typical flux-quantum 
densities n = BQic/e)" 1 w 2.5 -10 15 m~ 2 (in B — 10 T) 
and cyclotron scales = HeB /mc « 1 meV of electrons 
in fractional quantum Hall states. The excess density 
of Cooper pairs (beyond that of the nearest insulator) 
is controlled by the gate voltage in the superconduct- 
ing state, and can be brought near and below n$ to 
stabilize a fractional incompressible quantum liquid in a 
finite parameter range of size w$ surrounding the quan- 
tum critical point in Fig|2ta). Detecting fractional charge 
and statistics in the absence of magnetic fields, by exist- 
ing shot-noise or quantum interferometry methods from 
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FQHE experiments 3l|, |32j , would provide clear evidence 
of an established fractional TI in the quantum well. 

Without microscopic modeling and experimental data 
we cannot rule out a possibility that singlets would con- 
dense before triplets in a particular device. But even 
then, a further raise of the gate voltage would eventually 
condense triplets. Singlets cannot completely screen out 
the gate from triplets because the quantum well is very 
thin, and they repel each other stronger than they re- 
pel the triplets by the Pauli exclusion principle. Future 
experimental probes of topological spin dynamics may 
be able to reveal fractional 77 vortex liquids even if they 
coexisted with a singlet superconducting state of the <f> 
fields (which cannot screen spin). 

The striking features of these unconventional insula- 
tors are excitations with fractional charge, spin and ex- 
change statistics 33|. Their topological order exhibits 
chiral edge states and ground state degeneracy without 
symmetry breaking in non-simply connected spaces. The 
simplest fractional TIs arise from an extension of the 
Laughlin sequence to quantum spin-Hall states with con- 
served spin. They are characterized by two integers m± 
that satisfy ±m± > in the absence of U(l) magnetic 
fields. The charge jo and spin Jq density of the spin 
S z = ±1 Cooper pairs in Laughlin states are related to 
the SU(2) "magnetic" flux density $: 



Jo = 



* /J 1 



2ir \m + 



p — 



1 ; 1 

to + m_ 



The TR-symmetric demagnetized ground states also sat- 
isfy to + = — m_. Excitations can carry fractional charge 
SQ = n + /m + +n-/m- and spin SS Z = 0,4-/772+— n_/m_ 
specified by two additional integers n± , which also affect 
the exchange statistics. The transverse charge and spin 
conductivities of rj± in electric fields Ei = —eij(j)j are 
a% y = 3i/4>i = (ml 1 + mZ 1 )^ and a s xy = J(/ 4>i = 
(ml 1 — ml 1 )^ respectively. Hierarchical and non- 
Abelian spin-Hall states are also possible [U [2? 



The above type of fractionalization actually applies to 
a model TI whose gauge-theory description (|2J has a dif- 
ferent SU(2) gauge field than discussed so far, possibly 
realizable in ultracold atomic gases. The proper spin-like 
quantum number in a solid-state TI is the in-plane spin 
projected on the direction perpendicular to momentum, 
the eigenvalue of z(S x p) as evident from (JTJ) . If it were 
conserved, measuring it on a fractional TPs quasiparticle 
in the momentum p eigenstate would yield a fraction of 
ih/2. The spin quantum h would become fractionalized 
just like the charge quantum e. 

However, spin is not conserved and we must ask in 
what sense can the spin fractionalization occur. We 
can presently give only a partial but reassuring an- 
swer. There is no fundamental obstacle to fractional- 
izing charge because it is conserved and exhibits "cy- 
clotron" dynamics in the spin-orbit SU(2) gauge field 



as revealed by ((4]). Thus, TR-invariant fractional TIs 
can exist. Since combining an appropriate integer num- 
ber of fractional quasiparticles must reconstitute a single 
electron, the quasiparticles must inherit from electrons a 
spin-like degree of freedom. The TR-invariance requires 
that this degree of freedom transform like spin under TR 
and span two degenerate states. Its fractional quantiza- 
tion is guaranteed by the fundamental properties of vor- 
tex dynamics in incompressible quantum liquids, and its 
existence defines new types of topological order not yet 
encountered in FQHE systems. However, this "fractional 
spin" need not be a good quantum number. Depending 
on the nature of the SU(2) flux and perturbations such 
as disorder, the fractional quantum number generally has 
a combined spin and orbital character. 

In conclusion, we argued that fractional TIs of Cooper 
pairs can be obtained in TI quantum wells with induced 
pairing interaction. The realization of such states in het- 
erostructure devices rests upon the ability to achieve a 
combination of sufficiently large spin-orbit coupling, high 
mobility and low temperature, but the available materi- 
als and technologies offer very good prospects. 
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